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Abstract 

This paper describes salient features of new sets of parton distributions obtained 
by the CTEQ Collaboration^ based on a comprehensive QCD global analysis of all 
available data. The accuracy of the new data on deep inelastic scattering structure 
functions obtained by the very high statistics NMC and CCFR experiments provides 
unprecedented sensitivity to the flavor dependence of the sea-quark distributions. 
In addition to much better determination of the small x dependence of all parton 
distributions, we found: (i) the strange quark distribution is much softer than the 
non-strange sea quarks and rises above the latter at small-x; and (ii) the difference 
d — u changes sign as a function of x. A few alternative sets of viable distributions 
with conventional assumptions are also discussed. 
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Introduction In the current theoretical framework, high energy lepton-hadron and 
hadron-hadron interaction cross-sections a, both in standard model and in new physics 
processes, are related to calculable fundamental parton interaction cross-sections a by the 
QCD factorization theorems|I[] as a sum of integrals convoluting the latter with universal 
parton distribution functions. The parton distributions can, in principle, be determined 
from analyzing a set of standard experiments — deep inelastic scattering (DIS), lepton- 
pair production (LPP), high-p t direct-photon production, W- and Z-production, high-p t 
jet-production, heavy flavor production, ... etc. As both theory and experiments have 
matured and grown in complexity and as the scope of the standard experiments expands to 
various collider processes, the task of proper and systematic global QCD analyses requires 
coordinated efforts of both theorists and experimentalists familiar with details of both 
sides of such an analysis. For this reason, the CTEQ Collaboration has been developing 
the necessary tools for carrying out up-to-date global analyses based on previous work of 
Duke-Owens |2[] and Morfin-Tung||. 

Previous global analyses were based on data from the SLAC-MIT, EMC, CDHSW, 
and BCDMS deep inelastic scattering experiments; the E288, E605 lepton-pair production 
experiments; and the WA70 direct photon production experiment. [ffl Recently released 
NMC data[| on F 2 n /F 2 p , F 2 P - F™, and F%' d using a muon beam and CCFR data@ on 
F 2 F | using (anti-) neutrinos are expected to have a significant impact on QCD global 
analyses because of their extended kinematic coverage (particularly at small x), their 
high statistics and minimal systematic errors. The precision of the current generation 
of DIS experiments now far exceeds the size of next-to-leading order QCD contributions 
to these processes; thus they probe the full complexity of QCD mixing effects between 
quarks and gluons in a properly conducted QCD analysis. We find that these new accurate 
DIS data, supplemented by constraints imposed by data on lepton-pair and direct photon 
production, leads to substantially increased sensitivity to the flavor content of the sea- 
quark distributions as well as to the gluon distribution — features conventionally thought 
to be beyond the reach of analyses involving totally inclusive data. 

To demonstrate the need for new comprehensive global analyses and the potential 
for discovering new features of parton distributions in light of the new data, we show in 
table I values of x 2 P er degree of freedom obtained by comparison of old and new data 
with theoretical numbers obtained from parton distributions widely used in the current 
literature: KMRS0, MT-S1, MT-B20, and the recent MRS D0§ (which is based on 
earlier preliminary NMC and CCFR data). We note that all old sets do rather badly 
with the new data. Only the MT-B2 set is close to being acceptable, but still misses the 
CCFR data by 1.5 - 2 standard deviations. Even the MRS DO set (which used the new 
data as input) does poorly, with most of its high \ 2 a ls° coming from the CCFR data.0 

2 These large % 2 values are consistent with the original MRS analysis because the errors were enlarged 
in that analysis. 
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General Description of New Analysis For the current analysis, we use the following 

sets of data: BCDMS & Fg d , NMC Ff , & F 2 ' m /Ff ,0 CCFR F" Fe & Ff e DIS 
structure functions, E605 LPP da/dydr, and WA70, E706 & UA6 direct photon da/dydp t . 
With kinematic cuts (Q 2 > 4 GW 2 in DIS and LPP, W 2 > 12 GW 2 in DIS, and p t > 
4 GeV^ in direct photon production), the total number of data points used in most of our 
analysis is 917. We have tried higher Q 2 -cuts and found the results are not sensitive to the 
choice. For the very high-statistics DIS experiments, statistical errors have been reduced 
to such a low level that systematic errors often dominate the experimental uncertainties. 
The systematic errors clearly cannot be neglected as has been done in some past analyses. 
On the other hand, a complete treatment of these errors, including all relevant correlations 
(as is usually done in single-experiment analyses), would be (a) entirely impractical in 
the context of a large-scale global analysis, and (b) uncertain in statistical significance 
since such a diverse set of experimental data is involved. Thus, we adopt the common 
practice of combining the statistical and systematic errors in quadrature point-to-point. 
We have also tried to add the errors linearly, and found similar conclusions (with, of 
course, reduced absolute values of \ 2 )- The relative normalization of the various data sets 
is allowed to vary around 1.0 during the fitting process with an associated \ 2 included 
using the appropriate quoted experimental normalization errors. Heavy target corrections 
are applied to the neutrino structure functions to generate equivalent isoscalar structure 
functions for nucleons. The correction factor is taken from measured ratios of structure 
functions on light to heavy targets measured in muon scattering experiments. In the 
crucial low-x region, the default correction factor taken for this analysis is a conservative 
factor based on the assumption of saturation of the shadowing effect for A > 40. Since this 
correction is significant compared to current experimental errors, we have also investigated 
changes resulting from using a correction factor based on a fit of shadowing data off He, 
C and Ca with subsequent extrapolation to Fe. The differences will be noted later. 

We perform least \ 2 fits using primarily the MINUIT program. All results have been 
verified and investigated by a completely independent program package developed by 
Duke and Owens. We use the overall \ 2 as we H as individual x 2 ' s f°r each experiment 
as measures of the relative "goodness of fit" , but do not attach strict statistical significance 
to the absolute values of the x 2 ' s f° r reasons mentioned above. For quantitative next- 
to-leading order (NLO) QCD global analyses based on data of high accuracy, one must 
take into account a number of non-trivial theoretical and experimental considerations. 
These issues have been described in detail in several reviews, || , [|K| , |TT| hence will not be 
discussed here. 

3 Although the ratio measurement is not totally independent from the individual F% F & F% , it is 
in fact more significant than the latter since it involves considerably smaller systematic errors and a 
somewhat extended kinematic range; hence this data set is included upon the advice of some members 
of the NMC collaboration. 
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Recent experiments have called into question a number of traditional assumptions 
about the unknown input parton distributions (at some fixed scale Qq) such as SU(3) or 
SU(2) flavor symmetry for the sea-quarks. Since the diversity of data sets used in the 
analysis and the measurement accuracy have increased sufficiently to provide stringent 
constraints on the theoretical unknowns, we find it desirable to try fits which are free 
from most conventional prejudices and compare the results with those obtained with 
more traditional assumptions. This comparison turns out to be very revealing. 

For most of our analyses, the QCD evolution kernel, the effective coupling, and all 
hard scattering cross-sections are taken to be in the MS renormalization scheme; hence the 
extracted parton distributions are, by definition, in the MS scheme. For the convenience 
of certain applications, we also obtained parton distribution functions defined in the 
DIS scheme. These distributions are generated by independent fits to the same data.0 
Likewise, we generate representative leading-order fits for applications which use LO hard 
matrix elements. For concreteness, in the following description of our results, we shall 
focus on the MS distributions. 

Description of Results The most notable fact of this new global analysis is the ex- 
traordinary quantitative agreement of the NLO-QCD parton framework with the very 
high statistics DIS experiments over the entire kinematic range these experiments cover 
and the consistency of this framework with all available experiments in lepton pair and 
direct photon production as well. In the least constrained fits, the overall x 2 is typically 
around 860 for about 880 degrees of freedom (917 data points with 30 ~ 35 parton shape 
and relative experimental normalization fitting parameters) ; and the x 2 for the individual 
experiments are all comparable to the respective number of points — a remarkable result 
considering the diversity of processes, the kinematic range covered, and the accuracy of 
many of the experiments. Other fits with certain restrictions (such as SU(2) or SU(3) 
symmetric sea) have overall x 2 around 910 or higher. Although the relative values of x 2 
are considerably larger than the best fit, even these alternatives can represent reasonable 
fits when taken by themselves. The plots of these fits against data (not included in this 
letter) give the impression of good agreement in most cases and only very careful exam- 
ination will reveal the differences. The merit of these possible solutions can be better 
judged by examining details of the parton distributions obtained and of the comparison 
with the various data sets used, as will be discussed later. 

In general we find two common features of the new fits: 
(i) The small-x (~ 0.01 — 0.08) data from NMC and CCFR require substantially increased 
sea-quark distributions in this region. We found this also produces an indirect effect on 

This is a more reliable procedure then by transforming from the MS sets since the "NLO" term in 
the conversion formula can, under some circumstances, be comparable in size to the "LO" term, thus 
vitiating the perturbative nature of the transformation. S . 
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the shape of the gluon distribution such that the momentum fraction carried by the gluon 
at Q = 2 GeV is reduced from the previously accepted value ~ 0.45 to around 0.42. This 
fraction is partially recovered after evolution to larger Q . (Note that in NLO QCD the 
gluon momentum fraction is factorization scheme dependent (i.e. different in the MS and 
DIS schemes); and it is not exactly tied to the integral of the neutrino structure function 
F 2 as in leading order.) 

(ii) As a free-parameter in the global fit, the NLO Aq C d(5 flavors) is found to be typi- 
cally 150 to 160 MeV in these new fitsQ — consistent with typical values [Aqc<d(4/Z) = 
220 to 240 MeV ] found in direct d\n.F 2 / dhiQ 2 analyses [0] and in previous global 
analyses. [PI, Fits with A.Qc D (5fl) fixed at 220 MeV (corresponding to a higher value 
of a s favored by some LEP measurements) have substantially higher x 2 — by about 30. 

In Table II we present a set of representative fits with different characteristics and list 
their overall x 2 as we U as the x 2 P er data points for the various data sets. The best fit 
in the MS scheme, with the least constraints, is designated CTEQ1M.[] For this set, the 
input distribution functions at Qq (= 2 GeV) for all the flavors are taken to be of the 
form: f(x, Qq) = A x Al (1 — x) A2 (1 + A 3 x Ai ) . The coefficients Ai are subjected to minimal 
restrictions in addition to the quark-number and momentum sum rules. Although, it may 
appear that there are many parameters to be determined, it turns out that the wealth of 
data provides surprisingly stable results for the best fit irrespective of where the starting 
point is taken to be in the parameter space. A glance across the first row of Table II 
reveals that all data sets are uniformly well-fit[], indicating a good deal of consistency. 
The most conspicuous feature of the CTEQ1M distributions is the rather large flavor 
dependence of the sea quarks — the u, d, and s distributions all differ substantially in 
shape. See Fig.l. In addition, because of the significantly increased s(x) in the small-x 
region, the normalization of the strange quark, as measured by the ratio of momentum 
fractions (integrated over the full x-range (0,1)) k =2s/(u + d), is around 0.9 instead 
of the often quoted value of 0.5.0] Since these results appear to be unusual, we will 
discuss them in some detail after we describe the other fits in order to provide a basis for 
comparison. 

CTEQ1D and CTEQ1L are DIS and LO distributions with the same assumptions as 
CTEQ1M. The DIS distribution gives equally good fit to data as the MS one, just as one 



5 This result supersedes a preliminary CTEQ value reported at the 1992 International High Energy 
Conference at Dallas JTo| which is inaccurate due to a computer error. 

6 We adopt the following label scheme for our parton distributions: CTEQnSX, where n is a version 
number; S designates the factorization scheme - M for MS, D for DIS, and L for LO; X is absent for the best 
fits, otherwise it distinguishes the alternative sets within a given scheme. Parametrizations of these distri- 
butions can be obtained by sending a request to Botts@MSUPA (bitnet) or Botts@MSUPA.PA.MSU.EDU 
(internet) . 

7 The only apparently "high" x 2 associated with the NMC hydrogen date is due to differences of the 
NMC and BCDMS data around the x — 0.225 bin where the more abundant and more accurate BCDMS 
data dominate the fit. 
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expects - a consistent change in factorization scheme should not change the physics. The 
LO fit has a substantially higher y 2 (by ~ 50) and associated normalization correction 
factor for the LPP data, confirming the need of the NLO formalism for precision QCD 
analyses. CTEQ1MS is a set of distributions with a "singular" gluon behaving like x~ a5 at 
Q ; its x 2 , comparable to the best fits, indicates that current data allow a range of different 
extrapolations to smaller x beyond 0.01. And CTEQ1ML is a set obtained with A(5/Z) 
fixed at 220 MeV (the "LEP event shape value"). All these distributions give reasonable 
fits to data with different characteristics (and some with associated limitations, such as 
the LO one). We shall not discuss them further in this short report. 

The row labeled M2 represents a fit assuming a SU(2) flavor symmetric sea. The overall 
X 2 is higher than that of CTEQ1M by 50, all data sets are uniformly well-fit as for the 
latter. However, it is associated with a normalization of the strange quark corresponding 
to k — 1.41 - too high a value from the physics point of view. In order to investigate 
the sensitivity of the combined data set to the strange quark normalization and shape, 
several fits with k fixed at 0.5 (according to conventional wisdom) are tried. The last 
row of Table II, labeled K5, represents the best fit one can get with this restriction. The 
overall x 2 increased dramatically by 88. If one requires the shape of s(x, Qo) to be the 
same as the non-strange sea distributions, the x 2 a l so increases very substantially. Thus, 
contrary to conventional expectations, we found this data set (consisting of only totally 
inclusive quantities) is sensitive to k and to the shape of s(x,Q ) , and it strongly prefers 
K ~ 0.9 or higher. 

To trace the origin of this surprising result, we note the principal difference between the 
CTEQ1M and the K5-fits in Table II is the increased x 2 values of the CCFR F% and the 
NMC F£ data sets. In leading order QCD, F^ (converted from i 7 ^ 6 using a heavy-target 
correction factor) is given by the straight sum of all quark flavors whereas F^ corresponds 
to the same sum weighted by the squared charges. One can therefore get a direct handle on 
the strange quark distribution by examining the equality ^F^SFf = xs(x, Q) + 0(a s ), 
valid at small Q where the charm distribution is small. Fig. 2 shows a plot of the quantity 
on the left-hand side of this equation at Q 2 = 5 GeV 2 using data from NMC and CCFR, 
and compares it with the strange quark distributions of CTEQ1M, K5, and MRS DO at 
the same Q 2 . The size of the experimental errors are superimposed on the ^.F% d — 
curve. It is clear from this plot that k ~ 0.5 parton distribution sets are not compatible 
with the new data; and that both the size and the shape of s(x, Q) in the best fit are 
driven directly by data. How can this conclusion be reconciled with the conventional 
low value of K obtained from leading order parton model analysis of dimuon production 
data? |23 First, we note that our k is defined at a fixed Q 2 (=4 GeV 2 ) and is integrated 
over the full x-range (0, 1); whereas the k determined from the dimuon experiments use 
data from a wide range of Q 2 extending from less than 1 GeV 2 for the lowest x values 
to well over 100 GeV 2 at moderate values of x. The two quantities are not the same as 
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k is Q- dependent for a non-SU(3) symmetric sea. Secondly, from Fig.l, we see that for 
x > 0.1 our s(x,Q) actually becomes considerably smaller then u and d as one expected. 
If ft(x m i n ) is defined as the ratio of momentum fractions integrated from some x m { n to 1, 
then k(0.1) « 0.5; it increases to k(0.01) ~ 0.85 and eventually to m 0.9 in the theoretical 
limit x m i n — * 0. Thus, the comparison with the dimuon experiment is also dependent on 
the effective x-range and the extrapolations used in going beyond that range. Finally, as 



has been pointed out J|15|| existing LO dimuon analyses are likely to be unreliable at small 
x due to the hitherto neglected contribution from the gluon-fusion mechanism which is 
nominally next-to- leading order but physically significant. 

The question of SU(2) symmetry of the non-strange sea can be investigated by com- 
paring the two fits CTEQ1M (non-symmetric, k = 0.9) and M2 (symmetric, k = 1.4) in 
Table II. The increase in the overall x 2 of 52 mainly comes from the CCFR F 2 ^ and E605 
data sets. We know that the pN LPP cross-section should be sensitive to u because of 
the dominance of uu scattering.^ When restrictions are imposed on the flavor dependence 
of the sea, it is the interplay between the CCFR, NMC and E605 data sets which causes 
the x 2 to increase. While M2 has an (arguably) acceptable x 2 i h does have an abnor- 
mally large strange quark component. If one restricts the size of s(x,Qq), say k < 1.0, 
any solution with an SU(2) symmetric sea will have an unacceptable high x 2 ■ We also 
note, in the preferred set, CTEQ1M, the u and d distributions cross each other around 
x ~ 0.06. This behavior contrasts with that of MRS DO where (d — u) is chosen to be 
positive definite. We have found that the behavior shown is preferred for all the fits we 
have tried whenever the sign of this difference is left free.0 In this connection, we mention 
there is some tantalizing experimental evidence for this behavior of (d — u) from E772 at 



Fermilab.jnj We know of no reason why this quantity need be positive definite. 
Details on the above mentioned fits will be described in a full-length paper. 

Uncertainties and Challenges There are uncertainties associated with this analysis 
some of which can be addressed. Concerning the heavy target correction in the small- x 
region, we found the use of alternative shadowing correction schemes (cf. Sec. 2) leads 
to larger upward corrections to F% , hence an even bigger strange-sea. One may also 
consider the effect of shadowing inside deuterium on relating deuterium to nucleon struc- 
ture functions used in the analysis. Recent theoretical studies indicate that this effect is 



less than 2-3% at the smallest-x value covered here.|17| This is small compared to the 
uncertainties of the heavy-target correction. Detailed work can be done to quantify both 
these uncertainties. We also note, current analyses of DIS neutrino structure functions 
use the NLO QCD formalism for massless quarks, including the small but non-negligible 



5 The CTEQ analysis uses the full E605 data set on da/dydr (120 points); whereas, it appears from the 



published literature that, the MRS analyses use only the integrated da/dr data for consistency checks) 13 

9 This may be partially responsible for the very high % 2 values on the CCFR data using the MRS DO 
distributions (cf. Table I). 
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s — > c contribution. The effect of the charm quark mass is taken into account by applying 
a "charm-threshold correction" to the data, not at the theoretical end. This is not fully 
satisfactory because the applied "correction" is, at best, a leading order one. Since the 
complete NLO theory for this transition now exists, [IS] the correct procedure is to include 



this in the theoretical expression and to compare with the uncorrected physical structure 
functions. We do not, however, expect either of these possible improvements to affect 
the main features of our results since, as we have demonstrated, they are driven rather 
directly by the new data. The only alternatives to the large strange sea at small-x would 
be: (i) there are additional unspecified systematic errors associated with the two major 
experiments at small-x; and (ii) the theoretical corrections due to charm-mass effects are 
much larger than expected and turn out to account for the bulk of the observed difference 
1^2 d — 3F.2 d in place of xs(x, Q). The last possibility is under investigation. 

Hadron collider data on W-, Z-production, lepton pair production, direct-photon pro- 
duction, jet production, and heavy flavor production are providing useful tests of the 
conventional QCD calculations. The expected increase in integrated luminosity by an 
order of magnitude at the Tevatron will make these measurements important sources of 
quantitative information on parton distributions. The new parton distributions obtained 
in our global analysis lead to increased cross-section for all the above mentioned collider 
processes in the regions dominated by small-x partons (low mass lepton pair, low p t direct 
photon, and heavy flavor production). This is gratifying since earlier parton distributions 
tend to give smaller cross-sections on all these processes when compared to preliminary 
CDF data. Input from new results obtained in the current run at the Tevatron will help 
probe the very small-x region; and will measure, in addition, different combinations of par- 
ton distributions than currently available. Given the tight constraints on these functions 
already found in this analysis, the new input will certainly provide significant quantitative 
tests of the consistency of the QCD parton framework and lead to even better determi- 
nation of the parton distributions. We also, of course, eagerly anticipate results from the 



HERA ep collider |19| which will provide a wealth of information on structure functions 
and other observables sensitive to parton distributions in a widely expanded kinematic 
range. 
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Note Added in Proof: 

Since this paper was submitted, a new leading-order QCD analysis of dimuon events 
in neutrino scattering has been published by the CCFR group (PRL 70, 134 (1993)), 
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reaffirming a smaller strange sea than the non-strange one. The difference in the measured 
s(x,Q), compared to that obtained from our global analysis, lies mainly in the small x 
region. We have been informed by the CCFR collaboration that preliminary results from 
a next-to-leading order analysis of the same data do not lead to substantial changes in 
these results. If this measurement of the strange sea is correct, then the difference of the 
neutrino and muon F 2 structure functions shown in Fig. 2 cannot be fully attributed to 
the strange quark. This suggests either additional sources of theoretical uncertainties so 
far neglected in all global analyses, or an incompatibility of the F 2 measurements of the 
CCFR and NMC experiments within the stated errors in the very small x region. To 
see how changes in experimental systematic errors might affect the situation, we have 
performed new fits fixing k = 0.5 and increasing the experimental systematic error in the 
smallest-a; bins by a multiplicative factor. Reasonable fits with x 2 in-between the best 
CTEQ fits and the "k5 fit" (cf. Table II) are obtained with a factor of order 2 for either 
of the experiments. Aside from readjustments among the sea-quark flavors, the shapes of 
the other parton distributions remain quite stable. Details will be described elsewhere. 

It is important to emphasize that the question about the strange quark should not 
obscure the main results of this paper: the determination of a new generation of parton 
distributions based on recent highly precise experimental measurements and phenomeno- 
logical advances. The progress made is evident in the x 2 comparison between Tables I 
and II. The new degree of accuracy is also responsible for bringing the potential conflict 
between the new F 2 measurements and the dimuon results to light. To resolve the current 
puzzle we need to reevaluate all theoretical uncertainties; improve the treatment of sys- 
tematic errors in the global analysis; and, finally, include the dimuon data in the analysis. 
These efforts are underway. At the same time, we urge the key experiments to re-examine 
their data in the small- £ region to see if our observed discrepancy persists under closer 
scrutiny. 

We thank various members of the CCFR Collaboration, particularly Sanjib Mishra 
and Michael Shaevitz for in-depth discussions on the strange quark issue. 



9 





Total x 2 


ccfF 2 


ccfF 3 


nmcH 


nmcD 


nmcR 


bcdH 


bcdD 


E605 


Dir 7 


KMRS BO 


3954 


22.16 


3.34 


6.11 


8.46 


1.29 


1.25 


1.40 


1.20 


1.02 


MT SI 


3786 


21.56 


4.32 


4.43 


6.38 


1.81 


1.57 


1.44 


1.00 


1.77 


MT B2 


1415 


3.10 


3.14 


2.21 


1.41 


1.27 


0.76 


1.05 


0.90 


2.46 


MRS DO 


1643 


5.31 


3.71 


1.47 


1.24 


1.27 


1.17 


1.21 


1.16 


0.96 



Table I: x 2 of some widely used parton distributions applied to current experimental 
data. The first column gives the total x 2 f° r 917 data points (all experiments). The 
other columns give the x 2 P er data point for the individual data sets: CCFR F 2 and F 3 ; 
NMC H, D, and F?/F$; BCDMS H and D; £605; and the combined direct photon 
production data (WA70, £706 and UA6). 





Total x 2 


ccfF 2 


ccfF 3 


nmcH 


nmcD 


nmcR 


bcdH 


bcdD 


E605 


Dir 7 


# Pts. 


917 


77 


77 


83 


83 


89 


168 


156 


119 


47 


CTEQ1M 


860 


0.68 


0.71 


1.25 


0.90 


1.34 


0.72 


1.12 


0.87 


0.73 


CTEQ1D 


861 


0.74 


0.90 


1.26 


0.91 


1.30 


0.65 


1.07 


0.83 


0.70 


CTEQ1MS 


861 


0.71 


0.79 


1.26 


0.89 


1.25 


0.71 


1.10 


0.87 


0.74 


CTEQ1ML 


892 


0.78 


0.80 


1.24 


0.86 


1.28 


0.82 


1.18 


0.83 


0.75 


CTEQ1L 


914 


0.89 


0.86 


1.03 


1.07 


1.50 


0.68 


0.99 


1.15 


0.73 


M2 
k5 


912 
948 


0.81 
0.96 


0.96 
0.87 


1.28 
1.26 


0.87 
1.00 


1.30 
1.28 


0.76 
0.72 


1.11 
1.14 


0.99 
1.09 


0.89 
0.86 



Table II: x 2 of various CTEQ global fits to current data. Column labels are the same 
as in Table I. The number of fitting parameters (parton distribution shape and relative 
experimental normalization) is of the order of 30 to 35; resulting in about 880 degrees of 
freedom. 



10 



Figure Captions 



Fig. 1 : Sea-quark distributions. The flavor labels on the left are for CTEQ1M; those 
on the right for MRS_D0. 

Fig. 2 : Comparison of experimentally measured |-F^ d — 3i^ d with xs(x, Qq) from two 
of our fits, CTEQ1M and K5, and from MRS_D0. Errors of the two experiments are 
added in quadrature and superimposed on the relevant curve to indicate the size of the 
experimental uncertainty. 

References 

[1] J. Collins, D. Soper and G. Sterman in Perturbative QCD, Ed. A. Mueller, World 
Scientific (1989). 

D. Duke and J. Owens, Phys. Rev. D30, 49 (1984). 
J.G. Morfin and W.K. Tung, Z. Phys. C52 13, (1991). 

For a review and detailed references, see J.F. Owens and W.K. Tung, "Parton Dis- 
tributions of Hadrons" in Ann. Rev. Nucl. Sci. (1992). 

NMC Collaboration, P.Amaudruz et al, Preprint CERN-PPE/92-124. 

CCFR Collaboration, Nevis Preprints, S.R.Mishra et al, #1459; W.C.Leung et al., 
#1460; and P.Z.Quintas et al, #1461. 

J. Kwiecinski, A.D. Martin, R.G. Roberts and W.J. Stirling, Phys. Rev., D42, 3645 
(1990) 

A. Martin, R. Roberts, and J. Stirling, Durham and Rutherford Preprint 
(DTP/92/16, RAL-92-021) 

We used an emperical formula based on NMC measurements at small x and SLAC 
measurements at large x. (Private communication from A. Bruell (NMC). 

Wu-Ki Tung in Proceedings of 1992 International Conference on High Energy Physics, 
Dallas (to be published). 

See Wu-Ki Tung, in Proceedings of Workshop on Hadron Structure Functions and 
Parton Distributions, Ed. Geesaman et.ai, Singapore World Scientific, (1990). 



A. C. Benvenuti et.ai, Phys. Lett, 223B 490 (1989); P. Z. Quintas et.al. in |TJ and 



in 



P.N. Harriman, A.D. Martin, R.G. Roberts & W.J. Stirling, Phys. Rev. D42, 798 
(1990). 



11 



[14] See, for instance, H. Abramowicz et.ai, Z. Phys., C15, 19 (1982); K. Lang et.ai, Z. 
Phys. C33, 483 (1987); M. Shaevitz, in [0 

[15] M.G. Aivazis, F. Olness & Wu-Ki Tung, Phys. Rev. Lett. 65, 2339 (1990). 

[16] P.L. McGaughey et.ai, Phys. Rev. Lett. 69, 1726 (1992). 

[17] B. Badelek & Jan Kwiecinski, Nucl.Phys., B370, 278 (1992); V.R. Zoller Z.Phys., 
C54, 425 (1992); M. Strikman in |HJ. 

[18] See Ref.0; and J.J. van der Bij & G.J. van Oldenborgh, Z. Phys. C51 477 (1991). 

[19] Proceedings of the Workshop Physics at HERA, Ed. W. Buchmuller and G. Ingelman, 
DESY, 1991. 



12 



